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Abstract - The so-called grand canonical catastrophe of the density fluctuations in the ideal Bose 
gas is shown to be a particular instance of the much more general phenomenon of condensation 
of fluctuations, taking place in a large system, in or out of equilibrium, when a single degree 
of freedom makes a macroscopic contribution to the fluctuations of an extensive quantity. The 
pathological character of the “catastrophe” is demystified by emphasizing the connection between 
experimental conditions and statistical ensembles, as demonstrated by the recent realization of 
photon condensation under grand canonical conditions. 


Introduction. — The grand canonical behavior of the 
density fluctuations of an ideal Bose gas (IBG) is known 
to contradict the intuitive notion that fluctuations ought 
to be suppressed by lowering the temperature. This is 
commonly referred to as the grand canonical catastro¬ 
phe [TH3], alluding to a failure of the grand canonical en¬ 
semble (GCE) for bosons in the condensed phase. Fur¬ 
thermore, the successful realization of Bose-Einstein con¬ 
densation (BEG) of ultracold atoms, which requires the 
canonical or microcanonical framework, has contributed 
to consolidate the notion of naivety and inadequacy of the 
GGE. Yet, the GGE results are admittedly exact. Hence, 
the alleged flaws of the GGE do not arise from internal in¬ 
consistencies, but from the supposedly unphysical nature 
of the conditions required for the GCE to apply. How¬ 
ever, with the recent experimental realization of BEG in 
a gas of photons [3], complemented shortly after by the 
investigation of the condensate fluctuations [S], the out¬ 
look has changed, because condensation in this case has 
been achieved under grand canonical conditions. This is 
an important development with implications on the issue 
of statistical ensembles, since the grand canonical catas¬ 
trophe has been proven to be a real physical effect. 

In this paper we approach the fluctuation problem in the 
framework of large deviation theory and we show that the 
observation of the grand canonical catastrophe is of inter- 
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est in a context much wider than the physics of bosons, as 
an instance of condensation of fluctuations EHH]. Conden¬ 
sation of fluctuations is a phenomenon related to but dis¬ 
tinct from the usual condensation which appears after tak¬ 
ing a thermodynamic average. The latter, to be referred to 
as condensation on average, belongs to the realm of typical 
behavior, while the former normally is a rare event. These 
different instances of condensation may appear jointly or 
disjointly. Of particular interest is the case of non inter¬ 
acting systems in which the distinction between the two 
becomes extreme and, therefore, most clear since without 
interaction there cannot be condensation on average and 
yet fluctuations may condense [7]. In this respect, the 
IBG framework offers a rich and flexible scenario allowing 
to cover both cases: of condensation of fluctuations with 
and without condensation on average. This paper aims 
to show that when the two types of condensation occur 
simultaneously, the grand canonical catastrophe, rather 
than being a pathology, is the observable manifestation of 
condensation of fluctuations. 

In order to explain what condensation of fluctuations is 
about, it is convenient to lay down the basic probabilis¬ 
tic structure: Gonsider a generic system with microscopic 
states w in a macrostate described by the statistical en¬ 
semble P{u!\J,V), where J stands for a set of control pa¬ 
rameters and V for the size of the system. The setting is 
very general: w could be a single event in sample space 
as well as a trajectory, in which case V would involve the 
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trajectory’s time length. Here, since we have in mind an 
equilibrium application, V will be taken as the system’s 
volume. Referring to this ensemble as the prior, let M{uj) 
be an extensive observable which scales like V. Now, the 
probability that A4{uj) takes the value M is given by 

P{M\J,V)='^P{u}\J,V)5{M-M). (1) 

U) 

Though straightforward, this turns out to be a statement 
rich of consequences since a relation is established between 
the observation of a fluctuation in the prior and the impo¬ 
sition of a constraint on the same system. In fact, using 
terminology from statistical mechanics, P(M|J, H) plays 
also the role of the partition function in the constrained 
ensemble 

Pc(w|M, J, H) = J, V)5{M - M), (2) 

obtained by conditioning the prior to the event AI(w) = 
M. If AI(a;) obeys a large deviation principle, the above 
duality between fluctuations and constraints is reflected 
into the twofold role of the rate function 

I(m\J) = — lim — lnP(M| J, H), m = M/V (3) 
V->oo V 

as the entity controlling fluctuations in the prior and as 
the free energy density in the constrained ensemble [9]. 

The formal framework is completed by the important 
observation that the task of computing the rate function 
can be simplified resorting to a soft version of the con¬ 
straint (for details see [3] or [7]), namely by applying an 
exponential bias on the prior in place of the hard S func¬ 
tion constraint. This yields the biased ensemble 

Pb(cc|s, J, H) = (4) 

K{s,J,V) 

where if (s, J, V) = is the biased par¬ 

tition function. Continuing with the language of statis¬ 
tical mechanics, the constrained ensemble m describes 
the system enclosed by walls isolating with respect to Ai, 
while the biased ensemble describes the system in con¬ 
tact with an A4-reservoir. Then, the biased free energy 
Y{s, J) = limy_^oo y lnif(s, J, V) and /(m| J) are related 
by the Legendre transformation 

/(m|J) =s*m-r(s*, J) (5) 


The implication of this duality is that rare fluctuations, 
difficult to observe, can be made typical by the implemen¬ 
tation of constraints [10]. Conversely, if constraints are 
difficult or even impossible to realize in practice, in prin¬ 
ciple constrained system could be studied through the fluc¬ 
tuations of unconstrained ones. The key point, for what 
follows, is that fluctuations in the unconstrained system 
explore the phase diagram of the constrained one. 

The probabilistic nature of the above structure makes 
it of wide applicability. Among the many applications, as 
anticipated above, of particular interest is the apparently 
trivial case of an uncorrelated prior. In fact, the impo¬ 
sition of a constraint induces correlations which, in turn, 
may cause phase transitions to occur in the constrained 
system. Perhaps, the best known example of this type 
is the ferromagnetic transition due to the spherical con¬ 
straint imposed on the Gaussian model m- Then, as a 
consequence of duality, the transition singularity appears 
also in the behavior of the prior’s fluctuations while the 
prior’s typical behavior, due to the absence of correlations, 
is smooth and featureless by construction. In this case 
there will be condensation of fluctuations without con¬ 
densation on average. This phenomenon has been studied 
in different contexts such as information theory [m, fi- 
nance m and statistical mechanics EHH], encompassing 
both equilibrium and out of equilibrium behavior. 

Grand canonical vs mean canonical ensemble. — 

The rest of the paper is devoted to the study of the impli¬ 
cations of the general structure expounded above in the 
context of the IBG in equilibrium. Focus will be on the 
boson number fluctuations. In order to make the pre¬ 
sentation as simple as possible, a uniform system in a 
d-dimensional box of volume V with periodic boundary 
conditions will be considered. The microstates lo are the 
sets of occupation numbers ui = {np} of the single particle 
momentum eigenstates. The energy function is separable 
'H{uj) = the single particle dispersion rela¬ 

tion is of the form ep = ap“, where a is a proportionality 
constant. For instance, for photons a = c velocity of light 
and a = I, while for particles with mass m, a = l/(2m) 
and a = 2. 

The statistical ensemble is determined by the system’s 
preparation protocol. Thus, in the GCE which applies 
when the system is put in contact with a thermal and a 
particle reservoir, the probability of a microstate is given 

by 


where s* is the biasing field such that 

{M)s* = m (6) 

and where {■)s stands for the average in the biased ensem¬ 
ble. Summarizing, through J(m|J) a link is established 
between fluctuations in the prior and typical behavior ei¬ 
ther in the constrained Pc{uj\M, J,V) or in the biased en¬ 
semble Pb(w|s*, J, Id), which are equivalent in the sense of 
ensemble theory. 


Pge(c^|/3,M,H) = 


ZscW,^i.V) 




(7) 


where J\f{u!) = i® ^be number function. For large 

enough V, the equation of state takes the form [14] 


1 , X 

P= T7-, -- + A VW, 


( 8 ) 


VI-z 

where p is the average density, z = is the fugacity, A 
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is the thermal wavelength m, 

,, 1 r, 


(9) 


is the Bose function with v = d/a and T is the Euler 
gamma function. In the above formula the Bose function 
collects the contributions to the density from the excited 
states, while the F-dependent term is due to the average 
occupation of the ground state. In the i/ > 1 case, to 
which we will restrict from now on, gv{^) = where C 
is the Riemann zeta function. Hence, in the V ^ oo limit 
the isotherms become superiorly bounded by the critical 
value 

PcW) = A-V(l)- ( 10 ) 

This means that, by taking the thermodynamic limit and 
adopting the GCE protocol, it is not possible to prepare 
the system with an average density above pc- In other 
words, BEG does not take place in the GCE, as it is well 
known from the example of black body radiation, corre¬ 
sponding to the GCE with z = 1. Occurrence of BEG, 
as condensation on average, requires number conservation 
which can be implemented either rigidly or softly. In the 
first case, the number of particles is fixed and the ensemble 
statistics is canonical 


In the second case the number of particles is fixed on av¬ 
erage [14] by inverting the density-fugacity relation. This 
amounts to introduce the new ensemble defined by 


P„,c(w|/3,P,E) =Pge(a;|/3,zM/), (12) 

where the function z*{P, p, V) is the unique root of Eq. (|8]) 
given by (see Appendix 1) 

( ln5-i(A‘^p), for p<pc, 

\nz*{P,p,V)= { -AE-i/-, for p = PC, (13) 

I -^/[y{p- pc)], for p>pc, 

with A = l-A^r(iy)sm(7riy)/7r]y‘'. In the following, GCE 
will denote only the ensemble controlled by (/3, z, E), while 
the ensemble controlled by (/3,p, E) will be referred to as 
the mean canonical ensemble (MCE), adopting terminol¬ 
ogy from the spherical model literature [16]. Though for¬ 
mally similar, these ensembles become non equivalent m 
in the thermodynamic limit: Dividing the density axis into 
the normal phase below pc and the condensed phase above 
Pc, the GCE and MCE overlap in the first one but not in 
the second, which is accessible to the average density only 
in the MCE. 

It should also be noted that the origin of BEG in the 
MCE is in the mean-field treatment of the correlations 
generated by number conservation, with p* = In z* play¬ 
ing a role akin to that of the internal Weiss field in the 
mean-field theory of ferromagnetism. Once this is under¬ 
stood, the distinction between MCE and GCE becomes 


quite clear, since in the first ensemble the absence of in¬ 
teraction is formal, while the second one is genuinely non¬ 
interacting. This is the same distinction existing between 
the mean spherical model and the Gaussian model of mag¬ 
netic systems [TT][TB]. The similarities between the spher¬ 
ical model and the IBG have been pointed out by Kac 
and Thompson in the paper of Ref. m- On the exper¬ 
imental side, the control on the boson number is one of 
the key elements in the realization of BEC in the labora¬ 
tory. In experiments with ultracold atoms condensation is 
achieved in (micro)canonical conditions with control pa¬ 
rameters (/?, N, V) by enclosing a definite number of atoms 
in optical traps |3|. Instead, in the work of Klaers et al. |4], 
photons condensation has been obtained operating with 
the MCE protocol, that is by keeping fixed the average 
density p [T5] . 

Therefore, it is of interest to study the density fluc¬ 
tuations in the two ensembles, since on the basis of the 
considerations made in the Introduction, in the GCE (or 
in the MCE normal phase) condensation of fluctuations 
without condensation on average is expected, while in the 
MCE condensed phase the two types of condensation are 
expected to appear simultaneously. 

Density fluctuations. — As a consequence of the 
equivalence of the two ensembles in the normal phase, the 
study of fluctuations can be unified and carried out in the 
MCE framework for both of them. 

Using Eq. (HI, the probability to find the value N of 
M for a given p reads 


P^,iN\p,p,V) =J2 PM/3, ( 14 ) 

UJ 


We are interested in checking whether M obeys a large 
deviation principle and, if so, to find the rate function. 
This requires to extract the E-independent part from 
Imc = “'^In-Fmc, evaluated for large E. Carrying out 
the algebra expounded in Appendix 2 and denoting the 
fluctuating density by a; = N/V, for x < pc one obtains 


2mc(a;|/3,/0, E) 


z*{x) 1 

xln-^ + - \nz X 
z*{p) V 

1 ^^ z,M^*(p),v) 

V Z,4/d,z*(x),V)’ 


(15) 


with the function z* defined by Eq. m- Instead, for 
X > Pc one finds the linear behavior 


21mc(a;|/3,p, E) = {pc -x) In z*{p)+Imc{pc\/3,p,V). (16) 


Here, two comments are in order. The first one is about 
duality, which appears in a simplified form in Eq. (II3 
since Imc is related to the partition function in the same 
ensemble Z^c{P,x,V) = Zgd/d, z*{x),V). The reason is 
that by taking the MCE as prior the bias needed to render 
typical the fluctuation x amounts just to shift the fugacity 
from z*{p) to z*{x), without changing the form of the 
ensemble. Hence, with respect to number fluctuations, the 
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Fig. 1: Rescaled large deviation function vs rescaled flnctuating density in the MCE: normal phase (blue triangles) with 
— In 2 * = 0.4 and v = 3/2, condensed phase (black circles) with p > pc (Color on line). It is a straightforward computation to 
show that the temperature dependence is rescaled away by the A'^ factor. 


MCE is self-dual. The second one is that the linear branch 
of Eq. (fTBl) corresponds to condensation of fluctuations. In 
fact, exponentiating and comparing with the probability 
of the ground state occupation number 

Po(no|/3,z,E) = (l-2)e"'>l"^ (17) 

there follows that for x > pc and up to normalization one 
can write 

P„e(7V|/3, p, E) = Po(iV-iVc|/?, z*{p), V)P^,{Nc\P, p, E),| 

(18) 

which shows that macroscopic fluctuations above thresh¬ 
old can occur only through the zero momentum occupa¬ 
tion number, while the contribution of the excited states 
is locked to Nc = Vpc- 

The next step is to analyze separately the above result 
when p is in the normal or in the condensed phase. We 
shall see that the E —>• oo limit yields different behaviors 
in the two cases. 

Normal phase. As follows from Eq. (m, in the normal 
phase z*{p) < 1. Furthermore, z*{p) is independent of E. 
Thus, letting E oo, there exists the rate function 

Imcix\l3,p)= lim Ii„c(p|/3,P, E) 

V->oo 

+/3[/mc(/3,a:) - f,nc{f3,p)],x < 

(pc - x) Inz*(p)I^c{pc\P,z*{p)), X > PC, 

where the MCE free energy density is given by [14] 
fmci/3,p) = -f3-'^X~‘^g,,+iiz*{p)). As the plot displayed 


in Fig. [Tj (curve with triangles) shows, this is a convex 
function with an isolated minimum at the typical value p, 
where it vanishes, which is strictly convex for x < pc and 
linear for x > pc- 

Hence, in the normal phase of the MCE (or in the GCE) 
when X exceeds the critical threshold fluctuations do con¬ 
dense in absence of condensation on average. Though in¬ 
teresting, this fluctuation phenomenon is doomed to re¬ 
main an hardly observable event, since it is exponentially 
suppressed for large E. In fact, due to the existence of the 
isolated minimum at p 

lim Pmc(a:|/3, P, V) = S(x - p). (20) 

V—>-oo 

As we shall see, this situation is drastically changed in 
the condensed phase where condensation of fluctuations 
coexists with condensation on average. 

Condensed phase. When p is in the condensed phase, 
from Eq. m follows that z*{p) becomes E dependent 
with limy_,.oo z*{p) = 1. Hence, from Eqs. (ITSl) and (IT6)) 
follows 

Imc{x\/3,p) = 

a:lnz*(a:) +/?[/mc(^,a;) - fmc{P,Pc)],x < 

0 , x> PC ^ ’ 

which shows that the large deviation function vanishes 
identically for fluctuations above threshold [5] (see Fig.[T|). 
This means that when both p and x are in the condensed 
phase, fluctuations behave sub-exponentially with respect 
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to V and that in order to get the probability of fluctuations 
V dependent terms must be retained. Putting together the 
first line of Eq. (1211) and lowest order terms from Eq. dni) 
one has 


Imc(a;|/3,P, V) = 


Imc(xl^,p),x < PC 


,X> PC- 
( 22 ) 


Hence, due to the 1/E factor in the second line of the 
above equation, the probability of a fluctuation for x > 
pc becomes V independent and, in the E —>■ oo limit, 
instead of the 6 function of Eq. (EOl) now one obtains the 
distribution 


{ 0, for X < Pc, 

exp{- pip/; } (23) 

(P-Pcf ’ ^ > PC- 

which can be recognized as the Kac function of Ref. [T] and 
which is also the statement of condensation of fluctuations, 
as it can be checked by comparison with Eq. (ITT)) . This is 
the main result in the paper. 

With such a broad distribution the distinction between 
typical and rare events gets blurred and the width of fluc¬ 
tuations becomes macroscopic, since it goes like the den¬ 
sity of the condensate 

(24) 

Such a regime of strong fluctuations persists and is even 
enhanced by lowering the temperature, since pc vanishes 
as /3 —>■ oo. Therefore, these features, which constitute the 
grand canonical catastrophe, are nothing but the manifes¬ 
tation of condensation of fluctuations which, in the con¬ 
densed phase, is as typical as condensation on average. 

Conclusions. ~ The possibility to observe such a 
strong fluctuation regime depends only on the realiza¬ 
tion of the physical conditions specific of the MCE. Even 
though the experiment by Klaers et al. [J] is modeled by 
bosons trapped in an harmonic potential, yet the concep¬ 
tual structure expounded above still holds suggesting that 
the large fluctuations of the condensate reported in Ref. [S] 
are due to condensation of fluctuations. 

Summarizing, density fluctuations of the ideal Bose gas 
in thermal equilibrium have been analyzed in the normal 
phase (which is equivalent to GCE) and in the condensed 
phase of the MCE. In the first case condensation of fluc¬ 
tuations takes place as a rare event and in absence of con¬ 
densation on average. In the second case the two types of 
condensation occur simultaneously, both as typical events 
and producing the phenomenology of the grand canonical 
catastrophe, which, therefore, should not be regarded as a 
pathology, but just as an observable feature of the MCE. 
The natural question is whether the conditions respon¬ 
sible of these macroscopic fluctuations could be realized 
in atomic systems. Typically, atoms confined in optical 
traps are in microcanonical conditions and, therefore, the 
overall density does not fluctuate. However, one could 


think to have such a large system that a subsystem could 
be regarded in a GCE with energy and particle reservoirs 
due to the rest of the system. Then, if the chemical po¬ 
tential of the subsystem is controlled through the density 

I xt of the outer system, it was shown by Ziff et al. [I] 
at limv_),oo P{,x\l3, Pe^t,V) = 5{x - pext) for pext in the 
irmal as well as in the condensed phase. Hence, there 
would be no grand canonical catastrophe. However, it has 
been pointed out by Schmitt et al. [5] that this prediction 
is not in conflict with their findings, because it pertains 
to experimental conditions in which the subsystem is in 
diffusive contact with the environment, which is not the 
case for the photon experiment. The point is that the ex¬ 
periment of Schmitt et al. takes place in genuine MCE 
conditions, meaning by this that the control parameter is 
the average density p of the system itself, as opposed to 
Pext hr the conditions envisaged by Ziff et al. Though the 
distinction may seem rather subtle, yet it is enough to pro¬ 
duce different statistical ensembles as demonstrated by the 
drastically different fluctuation regimes in the condensed 
phase. 

Appendix 1. — Derivation of Eq. (TT^ : If p < pc the 
first term in the right hand side of Eq. ([8]) can be neglected, 
trivially yielding the E independent solution in the first 
line of Eq. m- Instead, if p > pc one has e = — ln 2 <C I. 
Therefore, using the expansions z ~ 1 — e and m 


9 u{z) = p^(l) -hr(l - ^ 


Eq. ([H]) can be rewritten as 


P- PC = 


Ee 




(25) 


(26) 


from which, it is straightforward to obtain the second and 
third line of Eq. (US, under the assumption v < 2, i.e. 
d < 2a, with r(l — v) = 7r/[r(zz) sin(7ri/)]. Notice that 
the condition on v is fulfilled for particles with mass and 
with d < 4. In the experimental conditions of Refs, 
photons do behave as particles with an effective positive 
mass. 


Appendix 2. ~ Using the integral representation of 
the Kronecker 6 

the probability (03 takes the form 
P„,c(iV|/3,P,E) = 

Jc 27ri z' 

(28) 

where C is a complex contour enclosing the origin and 

$(a;|/3,z,E) = -xlnz+ -^Zgc(/3, z, E), (29) 

with 

Zgc(/3, z, E) = ln(l - z) + A"V+i(2:)- (30) 
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Hence, 

2^mc(a;|/3,p,H) = lnPi„c(^|/3, P, H) 

= $(x|/3,z*(p),H)+/^(x|/?,H),(31) 

where 

K{x\P, H) = -1 ln£ ^ ^gy$(x|, 3 .z'.y)_ (33) 

The saddle point of the integrand function is determined 
by z\ V) = 0, which yields the equation 

x=^{.N')z', (33) 

whose solution z*{P,x,V), therefore, depends on x ac¬ 
cording to Eq. m- Hence, z*{x) is independent of V for 
X < Pc and H-dependent for x > pc- In the first case, 
from Eq. (1321) follows straightforwardly 

K{x\l3,V) = -^{x\(3,z*{x),V) + ^\tlz*{x). (34) 

In the second one, since z*{x) = ~ 1 with a 

very weak dependence on x, the integral is dominated by 
the region where the Bose function is well approximated 

by [IS] 

9 v+i{z)'Zi gv+i{l) + gi,{l)\iiz'. (35) 

As a check, one can verify that the saddle point equation 
obtained with the above substitution coincides, in the re¬ 
gion of interest, with Eq. Then, using the above form 
of gv+i{z') the integral can be carried out exactly, yielding 

if(x|/3,H) = -A-V+i(l), (36) 

which is independent of x. Inserting these results into 
Eq. §3), one finds 

I,nc(a;|/3,p,H) = 

<P{x\l3,z*{p),V) - ^{x\P,z*{x),V) + ^hYZ*{x), 

(37) 

for X < Pc and 

I„,c(x|/3, p, H) = Hx\/3, z*{p), V) - A-V+i(l), (38) 

for x > pc, which coincide with Eqs. (ITS]) and (flSl) . 
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